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SUMMARY 

The mass exchange between Hamilton Harbour and Lake Ontario, 
through the Burlington Canal, was considered important in 
estimating dissolved oxygen budget of the harbour. A lake- 
harbour exchange computational method was developed to 
utilize the current meter data (September 1-13, 1975) at two 
levels on a single tower. During the study period, the 
canal water was found to be quasi-isothermal; consequently 
unidirectional flow existed. 

It was estimated that an average of 2.14 x lo^ m^.d"^ 
(25 m'-s~M of harbour water flowed into the lake, while 
0.77 X 10^ m^.d"^ (8 m^.s~M of lake water flowed towards 
the harbour. This accounts for the total and net daily 
exchange of 1.04% and 0.50% of harbour volume with net 
exchange being towards the lake. As more water leaves than 
enters the harbour, the lake-harbour exchange is considered 
important for maintaining and even improving the existing 
harbour water quality. The harbour dilution factor was 
estimated as 0.0019 per day. 
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INTRODUCTION 

Hamilton Harbour is a natural harbour of fresh water in the 
north western corner of Lake Ontario. It is triangular 
shaped, 8x5 km, with a mean depth of 13 m (see Figure 1) 
and a water volume of 2.8 x 10® m^. The harbour is con- 
nected to the lake by the Burlington Canal which is approx- 
imately 820 m X 107 m wide x 9,5 m deep. Several creeks 
drain into the harbour with mean daily flow of 0.35 x 10^ m^ 
consequently, the natural throughput is approximately 0.12 
percent of the total harbour volume per day. The Hamilton 
and Burlington water pollution control plants discharge a 
total of 3.2 m^.s"^ of treated effluents. Together with the 
natural flow, this results in a total throughput of 0.22 
percent of the harbour volume per day. The major industries 
located on the south shores of the harbour use and recircu- 
late 27 m\s"* of harbour water (MOE, 1974). 

The harbour is hypertrophic, with mean chlorophyll a 20 
ug.£~^ and total phosphorus 40 to 200 yg.£~^. During the 
summer, the deeper waters of hypolimnion experience severe 
oxygen depletion. The lake-harbour exchange is important in 
determining a dissolved oxygen (DO) budget for the harbour. 
Based upon current measurements in the Burlington Canal, 
Palmer and Poulton (1976) stated that the mean exchange is 
of the order of 1% of the harbour volume per day. In the 
summer, volumetric exchange of 1% of the harbour volume per 
day accounts for 19% of the DO input to the harbour (Polak 
and Haffner, 1977) . 
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THEORY 

Freeman et al (1974) found the lake-harbour exchange is 
partly due to the Helmholtz mode, and van de Kreeke (1976) 
has shown that harbours with restricted entrances such as 
Hamilton Harbour have faster currents in the entrance due to 
the restriction. Helmholtz resonance is a balance between 
the kinetic energy of the flow through the canal and the 
potential energy due to the change in water level within the 
harbour. The kinetic energy of the flow through the canal 
causes a unidirectional open channel flow. Dick and Marsalek 
(1973) attributed the exchange to: 

(a) Unidirectional flow in the canal caused by the 
Helmholtz mode (difference in lake and water 
levels) . 

(b) Densimetric flow caused by the thermal strati- 
fication; warm harbour water flowing to the lake 
in the top layer and colder lake water flowing to 
the harbour in the bottom layer, underneath the 
warmer water. 

The unidirectional flow due to Halmholtz mode, persists 
while the densimetric flow occurs in the summer season only. 
The analysis by Dick and Marsalek (197 3) showed that the 
unidirectional flow is the major source of mass exchange 
through the Burlington Canal. 

As the exchange process is complex and variable, only 
estimates of the exchange volume are possible. Better 
methods can be formulated for a better data set consisting 
of current and temperature profiles in the vertical and 
horizontal plane at three cross-sections. In the absence of 
such data, an episode method is presented using current 
measurements from a single tower at two different depths. 
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Currents are known to be generally homogenous across the 
channel except for the wall boundary layers (about 2 m from 
the wall) . The data may, therefore, be transposed to the 
middle of the channel. The method outlined below applies to 
the isothermal case only {as determined by bathycasts and 
temperature measurements) when current speed profile with 
depth approaches parabolic distribution, expressed as; 

U^ = kY (1) 

where : 

U is the velocity along the channel axis at a depth 
Y from bottom. 

k is a constant 

Utilizing the average current speed at two depths (6.1 ra and 
7.5 m from bottom in 9.5 m of water), it can be shown that 
the mean current speed of the vertical cross-section is 80% 
of the average measured at the two levels. Let Ui be the 
mean current speed persisting towards the harbour for time 
ti and let the water particle travel an excursion distance 
Xi towards the harbour in time ti, such that: 

Xi = Uixti (2) 

Then the corresponding flow (towards the harbour) qi through 
the cross-section is: 

qi = 0.8 X Ui X A X ti (3) 

where A = cross-sectional area of the channel. 

The water movement may now reverse its direction due to 
Helmholtz mode. Let U2 be the mean current speed persisting 
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towards the lake for a time t2 and let the water particle 
travel an excursion distance X2 towards the lake in time tz, 
such that: 

X2 = U2 >< tz (4) 

Then the corresponding flow (towards the lake) qa through 
the cross-section is: 

q2 = 0.8 X U2 X A X t2 (5) 

The excursion distances, lying between 50 and 1100 m, were 
divided into 50 m classes upto 4 00 m (about half the channel 
length) and into 100 m classes for excursions greater than 
400 m. The mean of these classes E. were ranked with Ei 
being the maximum (Ei = 1050 m) . For each E. , the flows qi 
emd qj are grouped and summed, such that: 

For t = 1 to n 

H- = S qi (6) 

L. = I q2 C^) 

% 

T. = H . + L ( 8 ) 

where : 



H. is the total flow towards the harbour, 
t 

L. is the total flow towards the lake 

T. is the total exchange through the canal, 
t 



The above flows are then accumulated as: 



Hq. = Z H (9) 

^ = l 

i. 

Lq. = E L (10) 

^=l 

Tq^ = Hq^ + Lq^ (11) 
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A critical excursion distance E can be evaluated by the 

c 

dissolved solid (DS) budget as well as the plot of E. 
against T.. It will be shown that the flows, corresponding 
to excursion distance less than E , do not contribute to the 



c 



mass exchange. 



Therefore, for i = c 



1 



and E . = E , , , T « 

t c > (12) 

Qi = Hq,. = Hq^ 

Q2 = Lq. = Lq J 

X- c 

where Qj and Qz are total flows towards the harbour and 
lake respectively. 

Equations (9) , (10) and (11) yield accumulated or total 
flows when flows corresponding to excursion distances less 
than E. are neglected. 

It is apparent that the mass exchange through the canal 
results in dilution of conservative contaminants (e.g. 
dissolved solids) . Such a dilution factor may be computed 
as : 



D = (C;j - (C^,^ + C^))/(0.5 X (W1+W2) xp) (13) 

where : 

P = Dilution factor per day 

C, = Estimated mass of conservative contaminant 

n 

transferred from the harbour to the lake by 
the exchange process during p days. 



C„ = Estimated mass of conservative contaminant 
transferred from the lake to the harbour 
by the exchange process during p days. 
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C. = Estimated mass of conservative contaminant 
added to the harbour by the industrial and 
municipal discharges during p days. 

Wi = Estimated mass of conservative contaminant 
present in the harbour, initially. 

W2 = Estimated mass of conservative contaminant 
present in the harbour after p days. 

METHOD 

Two self-recording current instruments were operated at 6.1 m 
and 7.5m from bottom on a single tower in Burlington Canal 
(9.5 m total depth) from 14 August to 13 September, 1975. 
As the compasses of the Geodyne-type instruments have 
historically malfunctioned in the Burlington Canal, only 
Plassey current instruments (Model MO 21) were deployed. 
Prior to the installation, the instruments were calibrated 
and checked in the laboratory and set to record current 
directions, water temperature and current magnitudes (speed 
integrated over 10 min.) every 10 minutes. 

DATA ANALYSIS 

Water temperature and current speed data were pre-whitened 
(numerically smoothed) after Blackman and Tukey (1959; p. 29, 

39, 174) using binomial weights after Panofsky and Brier 
(1968; p. 150) as follows: 

C = Nl /(ml (N-m)I) (14) 

m 

where Cm is the binomial coefficient and these converge on 
the shape of the normal distribution curve as N increases. 

For the present data at 10 minute intervals, N = 4 results 
in the best smoothing. The binomal weights are then com- 
puted as: 
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m 
s = c / ): c (15) 

mm m 
o 



where S are the binomial weights for smoothing the data. 
The current directions were not smoothed as the water flow 
was restricted along the canal axis. The data were divided 
into monthly records for August and September. 

Water Temperature 

Figure 2 shows two bathythermograph casts in the Burlington 
Canal on August 15 and September 5, 197 5. Thermal stratifi- 
cation on August 15 is extensive while conditions appear to 
be approaching isothermal on September 5, 197 5. During 
August, thermal stratification existed in the canal creating 
a thermal lense. During thermal stratification, water 
movement data at four or more different depths are considered 
necessary to compute mass exchange as water may be moving in 
opposite direction at different levels. The available data 
at two levels (1.4 m apart vertically) were therefore not 
sufficient to compute mass exchange when stratification 
existed. The September data during nearly isothermal 
(unidirectional open-channel flow) conditions were processed 
for the computation of mass exchange as follows. 

The frequencies of occurrence for hourly temperature values 
were computed at both levels along with mean, standard 
deviation, maximum and minimum. These statistics are 
presented in Table 1. A simple correlation coefficient 
between the two temperature series was computed to be 0.81 
and X^ for the two series was computed to be 11.60 with 11 
degrees of freedom (see Table 1). As X^ = 19.68 
> X^ii = 11.60, it was concluded that temperatures at the 
two levels were not significantly different providing veri- 
fication of quasi-isothermal conditions. The computer 
plotting of currents at both levels (see Fig. 4) indicated 
that the flows were in approximate phase, consequently the 
quasi-isothermal conditions prevailed in the Burlington 
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Canal. Cross-correlations of water temperatures between two 
levels were calculated for up to 30 h lags and is presented 
in Table 2. The mean temperature at upper and lower levels 
was 20.5°C and 19.8°C respectively with standard deviations 



of 1.80C and a.S'^C. 



Water Movements 



Two dimensional frequency of occurrence of current magnitude 
and direction at both levels were computed and tabulated 
(see Appendix for details) . Table 3 presents a summary of 
statistical results. At both levels, the resultant currents 
during September were going towards the lake. The resultant 
speed at the upper level was 3.27 cm.s"^ while at the lower 
level it was 1.76 cm.s"^. The arithmatic average speed at 
both levels was approximately 8 cm.s"^ while the maximum 
speed recorded at both levels was just under 45 cm.s"^. 
The persistence factor at upper level was higher (0.40) than 
at the lower level (0.22). 

The smoothed current speeds were resolved along and across 
the canal and then averaged over an hour. The hourly values 
of speed and temperature were subjected to the spectral 
analysis (see Table 4) by standard numberical techniques and 
the Banning of the coefficients after Blackman and Tukey 
(1959; p. 34, 171). Cross-coherences between currents and 
water temperatures between the two levels were computed and 
significant coherences (95% confidence) are summarized in 
Table 5. 

Exchange Calculations 

As the currents going towards the channel wall do not 
contribute to the flow, all currents falling in the shaded 
area of Figure 3 were neglected. Water movements in the 
shaded area ranged from 3 ato 10% of the time. These few 
instances may be due to the instantaneous direction readings 
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recorded by the Plessey instruments. Normally, the water 
movements in the canal should be parallel to the axis of the 
channel. The rest of the currents at both levels (90-97%) 
were then resolved along the channel axis. 

The resolved speeds were then computer plotted for both 
levels (see example in Figure 4). Points to the left of the 
Zero Speed Line (ZSL) indicates flows towards the harbour 
while the points to the right of the ZSL show flow towards 
the lake. The date and time are also printed. When currents 
are mainly towards the walls, no points are printed during 
that interval. The points on the printout were then joined 
manually, the dotted lines indicate interpolation in the 
absence of one or more data points. Figure 4 shows that the 
two instruments recorded flows in approximate phase. As 
expected, the flows shift to and from the harbour each time 
the time history plot intersects the ZSL (see Figure 4). 
Between two such convective crossings, the flow persists in 
one direction. Now, if the average speed during such an 
interval is computed and a time interval noted, it is 
possible to calculate the excursion distance (Xi or X2 ) a 
particle may travel during that episode. The volume of such 
flow (qi or qa) can be calculated from equations (3) or (5) 
and are presented in Appendix. This process was repeated 
for the complete record of 13 days. Computations according 
to equations (6) to (11) were carried out and tabulated in 
Table 6. Thus, total flows toward the harbour and toward 
the lake were calculated for a period of 13 days for all 
excursion distances. Flow through the canal for different 
excursion distances is also presented in Figure 5 showing 
that small excursions (<325 m) are responsible for a maximum 
of 16% of the total exchange through the canal. However, 
the small excursions cause water to move within the confines 
of the canal length and never leave the canal. Therefore, 
these small excursions are responsible for exchange much 
less than the indicated 16% of the total. Thus, flows 
associated with small excursions (<325 m) are neglected in 
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the final estimates. A summary of exchange results is 
presented in Table 9. During the first IJ days ol Septembei 
1975, the average flow rate of water from the lake to the 
harbour was 0.8 x 10^ m'.d"^ or 9 m^.s"^ while the average 
flow rate of water from the harbour to the lake was 
2.1 X 10^ m^.d"^ or 25 ra^.s"^. This represented a total 
exchange rate of 2.9 x 10^ m\d"^ or 34 m^.s"^ (1% of the 
harbour volume per day) . The mass exchange for the harbour 
is schemetically presented in Figure 6. 

Dissolved Solids Budget 

The dissolved solids (DS) mass in the harbour were estimated 
for August 15 and September 16, 1975, utilizing a modified 
Boyce (1973) method (Polak and Haffner, 1977) and the DS 
vertical profile at a number of stations in the harbour (see 
Figure 1). Dissolved solids were comuted as 0.65 x conduc- 
tivity. The harbour was divided into 500 x 500 m cells each 
having constant depth. All cells were divided into one 
metre layers down to the bottom. The nearest DS profile to 
a cell was utilized to compute the DS content of the layer 
and the cell. The contents of all layers and cells were 
then summed for the total content of the harbour (Wj in 
equation (13)). Similarly, the content of the harbour after 
32 days was also estimated (W2 in equation (13)). 

The total DS added to the harbour during the 32 days (August 
15 to September 16) were estimated from industrial and 
municipal discharge data. Similarly, the DS addition (or 
reduction) to the harbour due to mass exchange was computed. 
The results of the above computations are presented in 
tables 7 and 8 indicating that the mass balance of DS can be 
approximately achieved if flows corresponding to excursions 
less than 325 m are neglected. 

Knowing the average values of DS in the harbour as well as 
in the lake (just outside the Burlington Canal on either 
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side) and the estimated flows in the canal in either direction, 
it is possible to calculate Cj^^ and C^ of equation (13). 
C. , Wi and W2 had previously been estimated for 32 (P) 
days. Using these values in equation (13), the dilution 
factor D was computed to be 0.0019 per day. 

DISCUSSION 

Water Temperature 

Chi-square testing for temperatures at two levels showed 
homogeniety. Table 2 lists cross correlations of tempera- 
tures between two levels for up to 30 h lags. At zero lag, 
the correlation coefficient was 0.81. Cross-coherence of 
temperatures between the two levels indicate significant 
coherences for 15.0, 17.1 and 40.0 to 125.0 h. The average 
temperatures at both levels were plotted (see Appendix for 
some examples) every two hours and these graphs confirmed 
the isothermal conditions in the canal. It may therefore be 
concluded that quasi-isothermal conditions prevailed in the 
Burlington Canal and consequently an unidirectional flow 
regime existed. 

Water Movements 

At both levels, the resulting currents moved towards the 
lake. The currents were negligible (< 0.3 cm.s"') for 20% 
of the recording period. Water movements were towards the 
lake (resultant direction) for 28% of the time while towards 
the harbour (opposite to the resultant direction) for 7% of 
the time indicating the net water movement towards the lake. 

Autospectra of Water Movements and Temperature 

Autospectral analysis of the currents (Table 4) indicates 
similar significant periods at both levels. The 12.0 h 
period may be related to the semidiurnal or tidal motions 
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and this periodicity occurs at both levels along the canal 
axis and at the upper level only across the canal axis. The 
12.1 h period for the temperature (composite value at upper 
and lower levels) indicates the semi-diurnal or tidal effects 
on water temperatures. A period of 10.9 h was observed in 
the temperatues at the upper level and this may be due to 
the lakewide effects. Kohli and Palmer (1973) found similar 
periodicities at Lakeview, Lake Ontario. The first mode of 
Lake Ontario oscillations (5.0 h) is observed in the Burlington 
Canal at both levels for currents as well as water temperatures. 
Rao and Schwab (1974) computed the first mode of oscillations 
for Lake Ontario as 5.11 h while Rockwell (1966) obtained a 
value of 4.91 h. The second mode (3.11 h) as computed by 
Rao and Schwab (1974) may be observed as 3.2 h in the currents 
and temperatures at the upper level (Table 4). The period 
of 2.4 h observed in the current along the canal axes at 
both levels and for temperatures at the upper level may be 
due to the third free oscillations of Lake Ontario (Rockwell, 
1966 and Rao and Schwab, 1974), Helmholtz mode for Hamilton 
Harbour (Freeman et al, 1974), or a resonant cooscillation 
involving both. 

The periodicities of the measured data relate to the theo- 
retical lake, harbour and Helmholtz periodicities, confirm- 
ing that the data set used for the mass exchange computation 
portrays the expected periodic motion. 

Mass Exchange 

There exists a critical excursion distance below which the 
water movements in the canal do not clear the confines of 
the channel. In order to evalute such a critical excursion 
distance, a DS budget for the harbour was computed. Table 7 
shows that during the 32 days (August 15 to September 16, 
1975), 12872 T of DS must leave the harbour through the 
Burlington Canal. Table 8 shows that the maximum (11625 T) 
amount of DS moves through the canal from the harbour to the 
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lake if the flows corresponding to excursions <325 ra are 
ignored. Thus, the critical excursion distance may be about 
325 m. This critical distance was also confirmed by Figure 5 
which showed the excursions <325 m may, at the most, contribute 
16% of the total flow. Based on a critical excursion 
distance of 3 25 m, average flow towards the harbour was 
9 m'.s~' (0.27% of harbour volume per day) while towards the 
lake, it was 25 m^.s" (0.77% of harbour volume per day). 
This resulted in a total mass exchange of 34 m\s~^ or 1.04% 
of harbour volume per day (see Table 9) . The net flow 

3 — 1 

towards the lake was, therefore, 16 m .s or 0.5% of harbour 
volume per day. The mass exchange resulted in a dilution 
factor of 0.0019 per day; consequently the water quality in 
the harbour may improve. 

During the period of no thermal stratification, (September 
26 to October 6, 1973), Palmer and Poulton (1976) estimated 
a net exchange of 1% of the harbour volume per day towards 
the lake and a total mass exchange of 3% of the harbour 
volume per day. These results are similar to the present 
study but larger in magnitudes. The larger flow estimates 
of Palmer and Poulton (1976) may be attributed to their 
computational method (based on the arithmetic mean values) . 
The present episode method yielded relatively better flow 
estimates as checked by the DS budget. 

CONCLUSIONS 

The mass exchange between Lake Ontario and Hamilton Harbour 
through Burlington Canal was found to be a complex process. 
Unidirectional flow and unstratified water in the canal were 
the basic assumptions of the present study which was based 
on data for September, 197 5. It is estimated that on the 
average 2.14 x 10^ m^.d"^ of harbour water flowed into the 
lake, while 0.77 x 10^ m^.d"' of lake water flowed into the 
harbour. This accounts for a total exchange of 1.04% of 
harbour volume per day with a net exchange of 0.50% of the 
harbour volume per day towards the lake . 
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The lake harbour exchange resulted in a dilution of 0.0019 
per day. This is very important for maintaining and even 
improving the existing water quality in the harbour because 
the harbour water, with its higher dissolved solids content, 
is discharged to the lake while better quality and oxy- 
genated lake water flows into the harbour. 

As the calculations for mass exchange were computed for 13 
days only and since they are based on the simplifying 
assumptions of unidirectional flow, they provide only 
partial information for numerical modelling of the harbour 
and accounting for the dissolved oxygen stocks. Refinements 
would require collecting data for more points and under 
stratified conditions and possibly improving the computation 
technique. 
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TABLE 1 : Water Temperature Frequency, September 1975 

Burlington Canal, Hamilton Harbour, Lake Ontario 



Temperature Range 


Upper Level 


Lower Level 


OC 




% 


% 


12.00 - 


12.99 


0.00 


0.14 


13.00 - 


13.99 


0.15 


0.98 


14.00 - 


14.99 


0.31 


1.12 


15.00 - 


15.99 


2.16 


3.77 


16.00 - 


16.99 


2.16 


7.40 


17.00 ~ 


17.99 


5.40 


11.31 


18.00 - 


18.99 


6.33 


8.66 


19.00 - 


19.99 


12.96 


12.15 


20.00 - 


20.99 


28.55 


15.36 


21.00 - 


21.99 


21.45 


20.67 


22.00 - 


22.99 


18.06 


14.66 


23.00 - 


23.99 


2.01 


3.49 


24.00 - 


100.00 


0.46 


0.28 


TOTAL 




100.00 


100.00 


MEAN °C 




20.48 


19.85 


STD DEV °C 




1.78 


2.28 


MAX °C 




24.48 


24.06 


MIN °C 




13.77 


12.87 


SERIES LENGTH 


660 


718 


Corelation 


Co-efficient = 


= 0.81 




xj^ = 11.60 







.'. Temperatures at two levels are not statistically different, 



TABLE 2 : Cross-correlations Between Water Tempera lines 
at Two Depths (Upper and Lower) 



Lag 
(h) 




I 
2 
3 

4 
5 
6 

7 
8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 



Upper Leads 


Lower Leads 




Lower 


Upper 




0.812 


0.812 




0.800 


0.807 




0.804 


0.800 




0.801 


0.809 




0.791 


0.806 




0.785 


0.783 




0.786 


0.774 




0.783 


0.772 




0.774 


0.764 




0.768 


0.762 




0.755 


0.765 




0.755 


0.755 




0.761 


0.747 




0.756 


0.757 




0.741 


0.763 




0.730 


0.759 




0.724 


0.742 




0.719 


0.731 




0.718 


0.726 




0.710 


0.722 




0.697 


0.724 • 




0.687 


0.717 




0.686 


0.706 




0.684 


0.702 




0.679 


0.707 




0.666 


0.708 




0.648 


0.701 




0.650 


0.696 




0.657 


0.692 




0.648 


0.695 




0.633 


0.702 





TABLE 3: Summary of Current Statistics in Burlington Canal, Hamilton 
Harbour, Lake Ontario, September 1975 



UPPER LEVEL LOWER LEVEL 



Resultant direction ccwiing 
from 0° as north 

Resultant speed (cm. s ) 

Average speed (cm. s ) 

Maximum speed (cm. s ) 

Persistence factor 

* 
Period of negligible speed 

(% of recording period) 

Percentage of time going 
in direction of resultant 

Percentage of time going in 
direction opposite to resultant 

Total number of readings 

Interval of readings (mins) 

Total depth of the Canal (m) 

Depth of the meter operation 
from bottom (m) 



252 



255 



3.27 


1.76 


8.23 


7.94 


44.94 


44.23 


0.40 


0.22 


20 


20 



37 



3925 
10 

9.5 

7.5 



19 



4309 

10 
9.5 

6.1 



* -1 

<0.3 cm. s 



TABLE 4: Summary of Major Spectral Periods (hours) 
Burlington Canal, Hamilton Harbour 
Lake Ontario 
September 1975 
80% Confidence Level 



Upper Level Lower Level 



* * * * * * 
Currents Along the 12.0, 5.0, 3.2, 2.4, 12.0, 5.0, 3.4, 2.4 

Canal Axis 2. 3 

* * 

Currents Across the 12,0, 5.0, 3.2 3.9 

Canal Axis 

* * * * * * 
Temperature-Grand Mean 10.9, 5.0, 3.2, 2.4, 13.3, 5.5, 3.3 

2.2 



* 

95% Confidence Level 



TABLE 5: Summary of Coherences (95% Confidence) 
Burlington Canal, September 1975 

Upper (llIA) and Lower (1110) Meter Locations 



Significant Periods (h) 



Speed along the Canal 

Speed across the Canal 

Lower - speed along Canal 
and water temp 

Lower - speed across Canal 
and water temp 

Temp at Upper and Lower 
locations 



3.2 

NONE 

12.0-10.0,6.0-4.6,3.9-3.4,2.9 

NONE 

120.0-40.0,17.1-15.0* 



90% Confidence 



Table 6: Flow through Burlington Canal, Hamilton Harbour. 





Mean 


Flow X 10 m 




Accumulated Flow x 


10^ m^ 


Average Accumulated 


Flow/day 




Excursion 
Travelled 














X 10^ 


m .d 






to harbour 


to lake 


Exchange 


to harbour 


to lake 


Exchange 


to harbour 


to lake 


Total 


i 


E. (m) 

X 


H. 

1 


L. 

1 


T. 

1 


Hq, 


Lq. 


Tq, 






Exchange 


1 


1050 





2358 


2358 





2358 


2358 


0.00 


0.18 


0.18 


2 


950 





1508 


1508 





3866 


3866 


0.00 


0.30 


0.30 


3 


850 


661 


1447 


2108 


661 


5313 


5974 


0.05 


0.41 


0.46 


4 


750 


1790 


3053 


4843 


2451 


8366 


10817 


0.19 


0.64 


0.83 


5 


650 


2224 


6387 


8611 


4675 


14753 


19428 


0.36 


1.14 


1.50 


6 


550 


864 


3412 


4276 


5539 


18165 


23704 


0.43 


1.40 


1.82 


7 


450 


2507 


6832 


9339 


8046 


24907 


33043 


0.62 


1.92 


2.54 


8 


375 


881 


1554 


2435 


8927 


26551 


35478 


0.69 


2.04 


2.73 


9 


325 


1035 


1330 


2365 


9962 


27881 


37843 


0.77 


2.15 


2.91 


10 


275 


1834 


1110 


2944 


11796 


28991 


40787 


0.91 


2.23 


3.14 


11 


225 


1608 


1695 


3303 


13404 


30686 


44090 


1.03 


2.36 


3.39 


12 


175 


1265 


983 


2248 


14669 


31669 


46338 


1.13 


2.44 


3.57 


13 


125 


836 


530 


1366 


15505 


32199 


47704 


1.19 


2.48 


3.67 


14 


75 


546 


218 


764 


16051 


32417 


48468 


1.24 


2.49 


3.73 



Table 7: Dissolved Solids (DS) Budget Hamilton Harbour, August- 
September 1975. 



r 



Date 



Total DS in Harbour 



15 August 75 

16 September 75 



75049 T 
73847 T 



From 15 August to 16 September (32 days), net 
reduction of DS in harbour = 1202 T 

Addition of DS to harbour by industry = 7424 T 
Addition of DS to harbour by municipal 
effluents = 4246T 



Total 



12872 T 



To achieve mass balance, 12872 T of DS must leave harbour through 
Burlington Canal and enter Lake Ontario/ during the 32 days. 



TABLE 8 : Net Transfer of Dissolved Solids (DS) to Lake 
August 15 - September 16, 1975. 



Ei 
(m) 



Net flow per day to lake 



1050 
950 
850 
750 
650 
550 
450 
375 
325 
275 
225 
175 
125 
75 



X 10 m 



0.18 
0.30 
0.36 
0.45 
0.78 
0.97 
1.30 
1.35 
1.38 
1.32 
1.33 
1.31 
1.29 
1.25 



Cumulative 
DS to lake in 32 days 

(T) 



1516 

2527 

3033 

3791 

6571 

8171 

10951 

11372 

11625 

11120 

11204 

11035 

10867 

10530 



Table 9: Summary of Mass Exchange Results Burlington Canal, 

Hamilton Harbour, Lake Ontario September 1-13, 1975, 



For 13 days and excusions 


> 


325 m 






Total flow into harbour 


= 


9.96 X 10^ m^ 






Total flow into lake 


= 


27.88 X 10^ m^ 






Total volume of Hamilton Harbour 


s 


2.80 X 10^ m^- 






Average daily flow into harbour 


^ 


0.77 X 10^ m^.d"^ 

3 -I 
8.87 m .3 








= 


0.27% of harbour volume 


per 


day 


Average daily flow into lake 


= 


2.14 X 10^ m^.d"^ 
24.87 m-^.s"-^ 








=: 


0.77% of harbour volume 


per 


day 


Total daily exchange 


^ 


2.91 X 10^ m-'.d""^ 
33.69 m^.s"-^ 




II 




= 


1.04% of harbour volume 


per 


day 



> 
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FIGURE I : HAMILTON HARBOUR AND BURLINGTON CANAL 
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FIGURE 2 ■ TEMPERATURE PROFILE BURLINGTON CANAL LAKE ONTAR 10 , 1975 
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FIGURE 3 : BURLINGTON CANAL 
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FIGURE 4= TIME HISTORY OF CURRENTS RESOLVED ALONG BURLINGTON CANAL, LAKE ONTARIO, 1975 
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FIGURE 5- MASS BALANCE OF HAMILTON HARBOUR, LAKE ONTARIO 1-13 SEPTEMBER, 1975 
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FIGURE 6: MASS EXCHANGE THROUGH BURLINGTON CANAL, LAKE ONTARIO, SEPTEMBER 1975 



APPENDIX 

Two-dimensional frequency of occurrence of current speed and 
direction for the two levels are presented in Tables 1.01 
and 1.02. 

Table 1.0 3 shows the detailed computational procedure for 
mass exchange. 

The temperature profiles (based on recorded values) every 
two hours are presented in Figure 1.01 for September 5, 
1975. 

Figures 1.02 and 1.03 display the dissolved solid profiles 
in the harbour on August 15 and September 16, 1975. The 
estimate of DS content of the harbour was based on these 
profiles. 
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TABLE 1.03: Detailed Computation of Mass Exchange between Hamilton Harbour and Lake Ontario, 1975 



Sept 


Start 




TO HARBOUR 










TO 


LAKE 








TO 




UPPE 


R 


LOVIER 




AVG 


UPPER 


LOWER 


AVG 


HARBOUR 


LAKE 












1 






















■^^ ( 








75 


Time 


i 


"l 


^1 


i 


"l 


''l 


i 


"l 


'^l 


i 


"2 1 


^=2 


i 


"2 


^^2 


i 


^2 


^2 


^-1 


^2 


1 


1.23 




















la 


^ 


648 


21 


5 


630 


19.5 


5.5 


639 




656 




4.53 


8 


3 


144 


6 


2 


72 


7.0 


2.5 


108 




















107 






6.13 




















7 


4 


168 











7.0 


4.0 


168 




171 




6.33 


6 


3 


108 


11 


5 


330 


8.5 


4.0 


219 




















208 






8.43 


7 


5 


210 


7 


5 


210 


7.0 


5.0 


210 




















214 






9.43 




















10 


4 


240 


8 


8 


384 


9.0 


6.0 


312 




3 30 




11.23 


8 
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240 


8 


5 


240 


8.0 


5.0 


240 




















245 






12.53 




















10 
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420 


10 


6 


360 


10.0 


6.5 


390 
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10 
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180 


9 
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162 
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16.03 
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TABLE 1.03: Cont'd 



Sept 


Start 








TO 


HARBOUR 














TO 


LAKE 








TO 


UPPER 


LOWER 


AVG 
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1 
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LAKE 
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7 


7 


294 


6 


4 


144 


7.5 


5.5 


219 




252 




19.13 




















11 


8 


528 


9 


8 


432 


10.0 


8.0 


480 




490 




23.13 


6 


9 


324 


7 


6 


252 


6.5 


7.5 


288 




















318 




4 


0.43 
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